Surface and bulk criticality in midpoint percolation 
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The concept of midpoint percolation has recently been applied to characterize the double perco- 
lation transitions in negatively curved structures. Regular d-dimensional hypercubic lattices are in 
the present work investigated using the same concept. Specifically, the site-percolation transitions at 
the critical thresholds are investigated for dimensions up to d = 10 by means of the Leath algorithm. 
It is shown that the explicit inclusion of the boundaries provides a straightforward way to obtain 
critical indices, both for the bulk and surface parts. At and above the critical dimension d = 6, it is 
found that the percolation cluster contains only a finite number of surface points in the infinite-size 
limit. This is in accordance with the expectation from studies of lattices with negative curvature. 
It is also found that the number of surface points, reached by the percolation cluster in the infinite 
limit, approaches 2d for large dimensions d. We also note that the size dependence in proliferation 
of percolating clusters for d > 7 can be obtained by solely counting surface points of the midpoint 
cluster. 

PACS numbers: 64.60. ah,05.70.Jk 



I. INTRODUCTION 



Percolation has been extensively studied over the past 
several decades and remains as an active field of re- 
search. In addition to its intrinsic scientific value and 
its role as one of the basic models of critical phenom- 
ena it has contributed to improving our general un- 
derstanding of statistical physics on various geometries 
(see, e.g., Ref. Q), as well as developing efficient nu- 
merical algorithms (for a general introduction to perco- 
lation, see Ref. Q). A percolation transition can be man- 
ifested in many different ways: a common quantity used 
in studies of the transition is the largest cluster size, but 
many other quantities also give a clear signature of the 
transition, including the cluster size distribution 0, H| 
and the ratio between the first and second largest clus- 
ter sizes @. At the same time, percolation is often a 
question of connectivity, so one obvious question is then 
how many connections can be randomly broken before 
the system fails to percolate from one side of the sys- 
tem to the other, or from its middle to the boundary. 
From this viewpoint, a percolation phenomenon requires 
a surface in order to be meaningful, just as water has to 
percolate all the way through the coffee layer in a cof- 
fee percolator. Thus, one might ask how the percolation 
transition is reflected on the actual surface. Specifically, 
the fraction of surface points belonging to the percolat- 
ing cluster, Coo, can be written as P s ~ (p — Pc) 13 " for 
d-dimensional systems where p is the occupation prob- 
ability and p c is its critical value. This formulation is 
parallel to the bulk criticality, i.e., P ~ (p — p c )P where 
P is the fraction of bulk points belonging to . The ex- 
ponent p s is known to be 4/9 from conformal invariancc 



for d = 2 [3], and the mean-field value, valid for d > 6, 
is f3 s — 3/2 Q. For d = 3, Monte Carlo simulations 
have estimated p s /v — 0.9753(3) where v describes the 
correlation length as £ ~ (p — p c )~ v [I]- 

The midpoint percolation turns out to be a useful 
concept to understand percolation transitions in curved 
structures (loj . This earlier work focused on lattices with 
negative Gaussian curvatures and it was found that the 
percolation for such lattices contains two critical thresh- 
olds: the first one at which the number of points reached 
on the surface from the midpoint becomes finite in the 
limit of the large system size, and the second one, where 
this number becomes a finite fraction of the surface 
points. The reason for the occurrence of two percola- 
tion thresholds appears to be intimately related to the 
size of the surface: when the surface-volume ratio is fi- 
nite, as for the negatively curved lattices, there appear 
two thresholds. These two separate thresholds coalesce 
into a single one if the surface- volume ratio vanishes. In 
a d-dimensional system with size N, this ratio is to lead- 
ing order N~ x / d . Thus from this viewpoint, a negatively 
curved structure with a constant surface- volume ratio is 
infinite-dimensional. 

Let b be the average number of surface points reached 
from the midpoint in a <i-dimensional hypercubic lattice 
having a linear size L. The quantity b has a size depen- 
dence of the form b ~ L K at criticality. The exponent 
k can be related to the bulk exponent j3 and the surface 
exponent /3 S as follows: 
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This follows since the probability that the percolation 
cluster at criticality contains the midpoint, in the large- L 
limit, is independent of how many surface points it con- 
tains. Thus, the fraction between the number of surface 
points and the total number of points of the percolating 
cluster becomes L^/L^ 1 ~ £-' 3 /' y L~^ v , where i"* 3 ^ 
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is proportional to P s and L is proportional to P (see 
Sec. lIII Al below). Consequently, the exact value f3 s = 4/9 
for d = 2 implies that k = 9/16 = 0.5625 0. A second 
case where the exponent b can be obtained analytically 
is d = 6, for which Eq. ([lj gives n = with the known 
values for the exponents This implies that only a 
finite number of points on the surface are reached from 
the midpoint at criticality. A third example is the Cay ley 
tree with coordination number z, which is a negatively 
curved lattice and hence can be regarded as correspond- 
ing to d — co. This structure has the lower and upper 
critical thresholds at p c \ = l/(z—l) and p c i = 1, respec- 
tively [Toj . and there exist indefinitely m any percolat- 
ing clusters between these two thresholds [11]. One can 
calculate the average number of surface points reached 
from the midpoint at the lower threshold p c \ which gives 
b = zj{z — 1). Since there is no size dependence, this 
means that k = 0. The implication is that one only 
reaches a finite number of surface points from the mid- 
point for all dimensions d > 6. 

One might also imagine a gradual reduction in the 
magnitude of the curvature until the structure becomes 
an ordinary d-dimensional lattice. Since both of the 
thresholds coalesce for a regular lattice, p c = p c i = p C 2, 
we are left with the following possibilities at criticality: b 
can be either a positive constant as it is for p c \ in case of 
the Cayley tree, or an increasing function of L as it is at 
Pci- In this work, we find that the divergent behavior of 
b indeed becomes weaker as d increases, and reaches the 
limit of b = const, at the upper critical dimension, d = 6 
in accordance with the expectation above. We also show 
that the surface observable, b, provides us with a direct 
way of quantifying the percolation transition in general 
dimensions. In Sec. UH we describe the Leath algorithm 
used throughout this work. In Sec. IIII1 we analyze the 
numerical results and discuss the critical behavior. Fi- 
nally, Sec. II VI gives a summary. 



II. LEATH ALGORITHM 

The measurement of b is well suited for the Leath al- 
gorithm 0, HJ , which is basically a burning algorithm to 
create a cluster starting from the midpoint. The collec- 
tion of burned sites at each realization will be called a 
midpoint cluster, and abbreviated as C rn . It suffices to 
count the number of surface points which are contained in 
C m . This makes the algorithm in the present case easier 
to use than in most earlier works since we need no ex- 
plicit treatment to exclude the boundary. Let us consider 
a d-dimensional hypercubic lattice having the length of 
each edge as L = 2™ with an integer n. For convenience, 
we will impose the periodic boundary condition in this 
lattice. Each point will be assigned a coordinate with 
d components. Fixing one of the components, one gets 
a (d — l)-dimensional surface across the system which 
serves as our effective boundary layer: we do not allow 
the midpoint cluster to grow beyond this layer. In this 



way, we effectively get an odd number of lattice sites in 
each direction which uniquely defines the midpoint. 

As for the actual implementation of the burning al- 
gorithm, one may choose between two options: depth- 
first and breadth-first methods. The former can be easily 
made using recursion as follows. 

1. Assign the "boundary" state to the points consti- 
tuting the boundary layers, and set all the other 
points as "unexamined" . 

2. Set the midpoint as "occupied" . 

3. For each of the neighboring unexamined points of 
this occupied one, 

(a) with probability p, mark the neighbor as oc- 
cupied and repeat Step 3 with respect to this 
newly occupied point. 

(b) With probability 1 — p, mark the neighbor as 
"stopped" . 

The midpoint cluster consists of the resulting occupied 
points and its growth will stop when it is completely sur- 
rounded by points marked as either stopped or boundary. 

In spite of the ease of implementation, this method 
requires that the program remember the state of every 
lattice point, which severely restricts accessible system 
sizes. For this reason, the other option, i.e., the breadth- 
first method is better suited for the present context. The 
program is then only required to remember the outmost 
shell of C m . This can be implemented in the following 
way. 

1. Prepare a queue and an array, both initialized as 
empty. The queue stores actively burning points 
with its possible directions to proceed, while the 
array stores points which cannot be burned again, 
i.e., either "stopped" or currently "active". 

2. Add the midpoint to the queue, with every direc- 
tion allowed, and record the midpoint as well as its 
state as active in the array. 

3. Retrieve an element from the queue. The queue 
becomes shortened by one in length on the retrieval. 
Also delete this active point from the array. 

4. For every possible neighbor from the retrieved 
point, make a search in the array unless it belongs 
to the boundary layers. 

(a) If it is absent, 

i. with probability p, add this to the array 
as active. Also add this to the queue with 
preventing it from propagating back. 

ii. With probability 1 — p, add this to the 
array as "stopped" . 

(b) If it is found as active, tell the corresponding 
element in the queue not to propagate toward 
this direction. 
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(c) If it is found as stopped, this neighbor is not 
penetrable. Do nothing. 

5. Go to Step 3. 

Note that we need not predefine the connection struc- 
ture since one can easily compute neighboring coordi- 
nates from a given point in regular structures. As rec- 
ommended in earlier works [12h14| ]. we have employed 



hashing |15| in Step 4, but not included the data block- 
ing method [ijj nor the generation of a random number 
from each site index [lj, [l6j]. We choose the burning 
probability p from the previously known site-percolation 
thresholds, and the values used are tabulated in Ta- 
ble. |H Since p c decreases with d, the actually gener- 
ated cluster is much smaller than the effective lattice 
size in a high dimension. We also note how efficiently the 
Leath algorithm performs with respect to memory usage. 
For example, our implementation with 2 200 megabytes 
memory can readily simulate a case where boundaries 
are away from the midpoint by 31 lattice spacings in 
a six-dimensional hypercube around the known thresh- 
old, p c » 0.109017 [14]. If we loaded every lattice point 
on memory, the memory requirement should roughly 
amount to 2 integers, a significant portion of which 
would be simply redundant since the burning probabil- 
ity is low. The depth-first method has been used in this 
work only for simulating small system sizes, as well as 
performance checks of the breadth-first method. 



III. RESULTS 

Using the above algorithms, we have generated more 
than 10 6 samples for each L and d obtainable within our 
resources. Our numerical results are presented in Figs. [1] 
and [5] and Table HI After recording the number of surface 
points reached by each midpoint cluster, we carry out the 
following analysis. 



A. Bulk criticality 

Let us consider a lattice at an occupation probability 
p, where one or more clusters are distributed over the sys- 
tem. Among them, our algorithm always picks only one 
cluster, C m , the midpoint cluster. When p reaches p Cl a 
percolation cluster Coo will appear among many clusters 
in the system. It may or may not contain the midpoint. 
In case that the midpoint happens to be contained in Coo > 
our algorithm shows that the generated cluster, C rn , actu- 
ally percolates the system, i.e., reaches the boundary. If 
Coo does not contain the midpoint, on the other hand, our 
algorithm picks up only a nonpercolating cluster with no 
surface points. In short, the frequency of touching the 
boundary in our algorithm simply means the probabil- 
ity for a percolating cluster Coo to contain the midpoint. 
Moreover, since the midpoint is one of the lattice points 
in the bulk, the probability that it is contained within 
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FIG. 1: (Color online) Data obtained by the Leath algorithm 
from d — 2 to 7. (a) The average number of surface points, 
b, reached by the midpoint cluster, (b) The average probabil- 
ity for the midpoint cluster to reach the boundary, which is 
identified with P, the density of the percolating cluster, (c) 
The average number of surface points contained in a perco- 



are shown in all the figures, 
symbol sizes. 



but usually comparable to the 



the percolating cluster Coo is the same as for any other 
bulk point. This probability is given by the density of lat- 
tice points in Coo and consequently has the size scaling 
P ~ L d f- d = L~ ,3 / v where df is the fractal dimension 
of the percolating cluster. This means that one can es- 
timate [3/p by counting how frequently C m with a given 
size L reaches the boundary since when it touches, it does 
belong to the percolating cluster with certainty and vice 
versa. We generally expect a size scaling of an observable 
Q at criticality as 



Q = AL (l + BL~ 



+ C, 



(2) 
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where A, B, and C are L-independent constants and 
A > characterizes the correction to the leading scal- 
ing behavior. In case of P, we know that T = —P/v. 
Furthermore, since P(L — > oo) in fact vanishes, we can 
safely set C as zero in this case. Therefore, one arrives 
at 

PLP' v w A + B'L- A (3) 

with B' = AB. Since plotting PL / V against L~ A should 
give a straight line for correct exponents, we can de- 
termine the exponents from the goodness of fit. For 
d = 2, for example, we get P/v = 0.1044(2) together 
with A = 1.4(2), when choosing the chi-square signifi- 
cance level as a — 10%, and this result agrees well with 
the exact value, p/v = 5/48 « 0.1042 (see Ref. [H and 
references therein). We have employed this method for 
lattices of dimension 2 < d < 7 and the results are pre- 
sented in Table U [Fig. Hp))]. Whenever the data fail to 
exceed the chi-square critical value (CV) for a = 10%, 
we remove the smallest system size and then repeat the 
procedure. As seen in Table U the resolution of our data 
tends to deteriorate somewhat with increasing dimension 
(see e.g., the values for d = 5). This is presumably caused 
by the fact that the correction term in Eq. © is not 
enough to absorb all the deviations from the leading be- 
havior when L is too small. Since the deviations increase 
with the dimension, we are restricted to using fewer sizes 
of L at a larger dimension. To some extent, the deviation 
is also caused by the fact that the number of available 
data sets decreases with large P/v, making the statistics 
worse. Finally we note that the deviation of P/v from the 
predicted values in case of the upper critical dimension, 
d — 6, could possibly also be attributed to the logarith- 
mic correction [la |. 

B. Surface criticality 

As for the bulk critical exponent /3, we use the sur- 
face data in order to obtain k. This time the observable 
is b, the number of surface points reached by the mid- 
point cluster C m [Fig. QJa)]. As discussed in Sec. HI the 
constant C is expected to be a nonvanishing constant 
for d > 6 since it corresponds to the number of surface 
points reached in the infinite-size limit. Thus, according 
to Eq. @, we assume the size scaling as 

{b- C)L- K w A + B'L- A , (4) 

which is of the same form as Eq. ([3]) apart from the addi- 
tional constant, C. As a practical data-fitting procedure, 
we find the smallest C that gives a sufficiently high CV 
to pass a = 10%. This procedure yields K = 0.563(1) for 
d = 2, which is entirely consistent with the exact result 
k = 9/16 = 0.5625. The result for 2 < d < 7 are given in 
Table HJ One way to obtain P s /v is then to use Eq. (TTJ) 
and the determined values for P/v and k. However, we 
will here use this connection as a consistency check as 



shown in the last column in Table Q] Instead, we use 
the alternative method of counting the average number 
of surface points reached by the midpoint cluster which 
actually do reach the boundary. This means that we only 
sample over the cases when C m does percolate [Fig.fljc)]. 
For d — 2, this gives P s /v = 0.33345(10), which is close 
to the exact value, 1/3. Table Q] shows p s /v estimated in 
this way for 2 < d < 7. To our knowledge, the values for 
four and five dimensions are reported for the first time 
in this work. It is notable that for d = 6 both P s /v and 
P/v deviate somewhat from the theoretical predictions 
but are still inside the estimated bounds. Nonetheless, 
the consistency check P/v + p s /v + n = d—1 is well born 
out for 2 < d < 6, verifying the internal consistency of 
our method and analysis (the last column in Table HI . We 
stress that the three exponents P/v, p s /v, and k in the 
present work are all obtained by just counting the num- 
ber of surface points reached by the midpoint percolation 
cluster at criticality. 

The case of d = 7 in Table U is of special interest be- 
cause P/v + Ps/v + k = 5.1(2) < d — 1 = 6. This means 
that the relation L K /L d ~ l ~ L" /3 / l/ L~ /3s / l/ breaks down 
in this case. The reason is that the percolating cluster 
at criticality is no longer unique for d > 7 (27j. This is 
usually referred to as the breakdown of the hyperscal- 
ing relation ||. Suppose that the number of percola- 
tion cluster scales as L x , then the relation changes to 
L K /L d - X - L-x L -fi/y L -fis/y so that the consistency re- 
lation becomes P/v + P s /v + k = d — 1 — x. The reason 
is that the chance that the midpoint belongs to one of 
the percolating clusters is still given by L df ~ d , but now 
L—0/ v gives the probability that a lattice point belongs 
to any percolating cluster. Consequently, the chance for 
the midpoint to be contained in one of percolating clus- 
ters is L~ x L~fil l/ from which P/v + P s /v + n = d— 1 — x 
follows. The number of percolating clusters is expected 
to scale as L d ~ 6 for d > 6 (27|. The consistency relation 
for d > 6 then becomes P/v + P s /v + n = 5. As seen in 
Table U this relation is born out by our results for d = 7. 
We also observe that, provided that the mean-field values 
hold for the critical indices P/v and P a /v for d > 6, and 
in addition the critical index k — for d > 6 (as implied 
by our results), then it follows from the consistency re- 
lation that x — d — 6. Conversely, if we take the growth 
of the percolation clusters, L d ~ 6 , for granted, then we 
obtain k — from the consistency relation. This also 
implies that the fractal dimension of a single percolating 
cluster becomes df = d — x — P/v = 4 for d > 6 [27|- It is 
interesting to note that the exponent x, describing how 
the number of percolating clusters grows with a lattice 
size at criticality for d > 6 can be obtained by just count- 
ing the number of surface points reached by the midpoint 
cluster. 

Next we assume that the inferred result n = for d > 6 
is correct. Then Eq. ((4]) reduces to 

b^b OQ +B'L- A , (5) 

where b^ = A + C. This says that b(L —> oo) approaches 
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TABLE I: Occupation probabilities used in this work and values of exponents obtained in this work in comparison to previously 
known results. For d = 2 and d > 6, the known values of ft/v and /3 s /v are exact. The last column is for a consistency check 
of our analysis where the sum should be d — 1 according to Eq. UJ). The discrepancy at d — 7 signals the proliferation of 
percolating clusters, as explained in the text. 



d 


p c (known) 


p/v (X) 


P/v (known) 


K(Y) 


(Z) 


Pa/v (known) 


X+Y+Z 


2 


0.5927460 a 


0.1044(2) 


5/48° 


0.563(1) 


0.33345(10) 


1 /3 m 


1.001(1) 


3 


0.3116081 b 


0.478(2) 


0.474(6)' 
0.477(4) s 
4770('2~) h 
0.481(l) i 

0.4774(l) j 


0.554(7) 


0.974(2) 


0.970(6) f 

0.975(4)* 
0.9754(4) n 
0.9753(3) j 


2.01(1) 


4 


0.196889 c 


0.945(5) 


0.9528(14) b 
0.953(7) c 


0.410(1) 


1.64(2) 




3.00(2) 


5 


0.1407966 d 


1.5(1) 


1.462(16) k 
1.46(l) c 


0.11(7) 


2.408(5) 




4.0(2) 


6 


0.109017 d 


1.9(1) 


2 1 


0.06(2) 


2.8(3) 


3° 


4.8(4) 


7 


0.0889511 d 


2.0(1) 


2 1 


0.06(8) 


3.08(1) 


3° 


5.1(2) 



a Reference 
b Reference 
c Reference 
d Reference 
c Reference 
f Reference 
g Reference 
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n Reference [2 
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a constant denoted as boo. We have estimated this con- 
stant up to d = 10, using the critical probabilities re- 
ported in Ref. [14| [Fig.^a)]. The correction exponent A 
shows a tendency to increase as d grows so that the con- 
vergence to &oo becomes more rapid. Figure[2jb) suggests 
that the limiting value of b^ for large d could possibly 
be 2d. In order to investigate this further, we note that 
b'(L) = p c b(L — 1) corresponds to the number of points 
reached on the surface provided that percolation along 
the surface is prohibited. In Fig. H^c), we have plotted 
b'^ = b'(L — > oo) against l/(d — 6), where we assume 
that d — 6 is a fundamental parameter in the problem. 
Linear extrapolation suggests that the limiting value for 
d — > oo is close to b'^ = 1. Since p c approaches l/(2d— 1) 
for large d, this result also implies that b^ w 2d, even 
though the precision is not sufficient to make any firm 
conclusion. Note that 2d is just the number of faces in 
a d-dimensional hypercube, so that one can say that C m 
reaches every face of the surface in one point on the aver- 
age. However, any individual realization of the midpoint 
percolation cluster will of course reach a variety of points 
on each face and some not at all. 



In a previous study, it was shown that the cluster grows 
on the average by one site per step for large enough di- 
mensions, suggesting a connection to the self-avoiding 
random walk (SAW) [14J : A cluster that grows with pre- 
cisely one site per step traces out a SAW. From this 
perspective, it is interesting to note that our quantity 
b^ corresponds to the average number of points that a 
SAW walker starting from the midpoint reaches on the 
surface. Since a SAW walker always has a finite chance 
of getting stuck in any dimension d < oo, this means 
that b'ac < 1 for SAW, whereas we have found b'^ > 1 
for the midpoint percolation cluster at criticality. How- 
ever, the chance of getting stuck vanishes for the SAW 
walker in the limit of d — > oo, which means that b'^ = 1 
in this limit. This agrees with our corresponding result 
for percolation, as well as with the result in Ref. [l4j |. 
suggesting some similar feature between an unhindered 
SAW and the midpoint cluster at criticality in the limit 
of large d. However, the cluster created by SAW and 
the midpoint percolation cluster at criticality have quite 
different structures in all dimensions including the limit 
of d — > oo since the SAW cluster has fractal dimension 
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df = 2 for high c? whereas that of the midpoint perco- 
lation cluster is df — 4 as verified in the present paper. 
One may note that if one only takes the backbone of the 
percolating cluster, the backbone indeed has a fractal di- 
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IV. SUMMARY 
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FIG. 2: (Color online) (a) Measurements of b/d for dimensions 
from d — 7 to 10, where b is fitted by Eq. ([5]). The correction 
exponents are estimated as A(d = 7) = 0.5(1), A(d = 8) = 
0.8(2), A(d = 9) = 0.7(4), and A(d = 10) > 1. (b) The 
limiting number of surface points, b(L —¥ 00), divided by d. 
(c) The number of surface points reached, when percolation 
along the surface is prohibited, as a function of l/(d — 6). 
Linear extrapolation suggests that 6^(d — > 00) ~ 1. 



We studied the percolation transitions in dimensions 
from d — 2 to 10 by the Leath algorithm. In particular, 
it was shown that one can obtain various critical proper- 
ties by just counting the number of points reached on the 
surface from the midpoint. To this end, we checked that 
our midpoint percolation yielded consistent results with 
known ones, and estimated exponents characterizing the 
surface criticality as well as the bulk one. We found that 
the divergent behavior of b becomes weaker in higher di- 
mensions as anticipated from the Cayley tree, so that it 
scales as L 9 / 16 for d = 2 but converges to a constant in 
the infinite-size limit for d > 6. We also confirmed that 
the percolation cluster ceases to be unique for d > 6. In 
addition, our results suggests that the number of surface 
points reached approaches the value of 2d in the limit 
of large d. Provided that percolation along the surface 
excluded, this corresponds to the simple result that pre- 
cisely one surface point is reached at criticality. 



Acknowledgments 

S.K.B. and P.M. acknowledge the support from the 
Swedish Research Council with the Grant No. 621- 
2002-4135, and B.J.K. was supported by the Korea Re- 
search Foundation Grant funded by the Korean Govern- 
ment(MEST) with the Grant No. KRF-2008-005-J00703. 
This research was conducted using the resources of High 
Performance Computing Center North (HPC2N). 



[1] K. Christensen and N. R. Moloney, Complexity and Crit- 
icality (Imperial College Press, London, 2005). 

[2] S. N. Dorogovtsev and A. V. Goltsev, Rev. Mod. Phys. 
80, 1275 (2008). 

[3] D. Stauffer and A. Aharony, Introduction to Percolation 
Theory (Taylor & Francis, London, 2003), 2nd ed. 

[4] P. L. Leath, Phys. Rev. B 14, 5046 (1976). 

[5] P. L. Leath and G. R. Reich, J. Phys. C 11, 4017 (1978). 

[6] C. R. da Silva, M. L. Lyra, and G. M. Viswanathan, 
Phys. Rev. E 66, 056107 (2002). 

[7] C. Vanderzande and A. L. Stella, J. Phys. A 20, 3001 
(1987); C. Vanderzande, ibid. 21, 833 (1988); S. Roux 
and A. Hansen, Phys. Rev. B 38, 5170 (1988). 

[8] J. P. Carton, J. Phys. (France) Lett. 41, L175 (1980). 

[9] Y. Deng and H. W. J. Blote, Phys. Rev. E 72, 016126 
(2005). 



[10] S. K. Baek, P. Minnhagen, and B. J. Kim, Phys. Rev. E 

79, 011124 (2009). 
[11] I. Benjamini, R. Lyons, Y. Peres, and O. Schramm, Ann. 

Probab. 27, 1347 (1999). 
[12] C. D. Lorenz and R. M. Ziff, Phys. Rev. E 57, 230 (1998). 
[13] G. Paul, R. M. Ziff, and H. E. Stanley, Phys. Rev. E 64, 

026115 (2001). 
[14] P. Grassberger, Phys. Rev. E 67, 036101 (2003). 
[15] D. E. Knuth, The Art of Computer Programming, vol. 3 

(Addison- Wesley, Upper Saddle River, NJ, 1998), 2nd ed. 
[16] H. Vollmayr, J. Stat. Phys. 74, 919 (1994). 
[17] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 
[18] J. W. Essam, D. S. Gaunt, and A. J. Guttmann, J. Phys. 

A 11, 1983 (1978). 
[19] P. Grassberger, J. Phys. A 25, 5867 (1992). 
[20] Y. Deng and H. W. J. Blote, Phys. Rev. E 69, 066129 



7 



(2004) . 

[21] Y. Deng and H. W. J. Blote, Phys. Rev. E 71, 016117 

(2005) . 

[22] R. M. Ziff, Phys. Rev. Lett. 69, 2670 (1992). 

[23] H. G. Ballesteros, L. A. Fernandez, V. Martfn-Mayor, 

A. M. Sudupe, G. Parisi, and J. J. Ruiz-Lorenzo, Phys. 

Lett. B 400, 346 (1997). 
[24] H. G. Ballesteros, L. A. Fernandez, V. Martfn-Mayor, 

A. M. Sudupe, G. Parisi, and J. J. Ruiz-Lorenzo, J. Phys. 

A 32, 1 (1999). 
[25] J. Adler, Y. Meir, A. Aharony, and A. B. Harris, Phys. 



Rev. B 41, 9183 (1990). 
[26] A. B. Harris, T. C. Lubensky, W. K. Holcomb, and 

C. Dasgupta, Phys. Rev. Lett. 35, 327 (1975). 
[27] L. de Arcangelis, J. Phys. A 20, 3057 (1987). 
[28] B. Sapoval, M. Rosso, and J. F. Gouyet, J. Phys. (France) 

Lett. 46, L149 (1985). 
[29] We previously conjectured in Ref. [l(| that k — 4/7 « 

0.5714 for d — 2 by relating this exponent to the fractal 

dimension of the hull 28]. In the present work, we can 

distinguish between 4/7 and 9/16 and thus verify that 

the latter is correct. 



